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Abstract 



In the internal shock model of gamma ray bursts ultrahigh energy muons, pions, neutrons 
and kaons are likely to be produced in the interactions of shock accelerated relativistic protons 
with low energy photons (KeV-MeV). These particles subsequently decay to high energy neu- 
fvq ■ trinos/antineutrinos and other secondaries. In the high internal magnetic fields of gamma ray 

bursts, the ultrahigh energy charged particles (/i+, 7r+, K^) lose energy significantly due to 
synchrotron radiations before decaying into secondary high energy neutrinos and antineutrinos. 
The relativistic neutrons decay to high energy antineutrinos, protons and electrons. We have 
calculated the total neutrino flux (neutrino and antineutrino) considering the decay channels 
(~| , of ultrahigh energy muons, pions, neutrons and kaons. We have shown that the total neutrino 

O4 flux generated in neutron decay can be higher than that produced in /z"*" and tt'^ decay. The 

Q . charged kaons being heavier than pions, lose energy slowly and their secondary total neutrino 

$-H ■ flux is more than that from muons and pions at very high energy. Our detailed calculations on 

C/5 , secondary particle production in pj interactions give the total neutrino fluxes and their flavour 

ratios expected on earth. Depending on the values of the parameters (luminosity, Lorentz fac- 
tor, variability time, spectral indices and break energy in the photon spectrum) of a gamma ray 
CO ■ burst the contributions to the total neutrino flux from the decay of different particles (muon, 

^ , pion, neutron and kaon) may vary and they would also be reflected on the neutrino flavour 

^^ ■ ratios. 

.^ ■ Keywords: Gamma Ray Bursts, Neutrinos 

5 ■ 1 Introduction 

Extensive airshower arrays [1, 2, 3, 4, 5] have detected a large number of ultrahigh energy cosmic ray 
events. The origin of these cosmic rays are yet to be identified. They may come from supernova 
r^ . remnants (SNRs), active galactic nuclei (AGN), gamma ray bursts (GRBs) or some unknown 

C^ I sources. Some of the ultrahigh energy cosmic rays (protons and nuclei) may interact with low energy 

radiations and matter inside the source producing secondary high energy photons and neutrinos. 
Protons can be shock accelerated to ~ 10^^ eV inside GRBs by Fermi mechanism. Shock accelerated 
protons may loose energy by synchrotron cooling and p'j interaction inside the source depending 
on the magnetic field and low energy photon density respectively. The high energy neutrino flux 
produced in p^ interactions have been calculated in detail in many earlier papers [6, 7, 8, 9, 10, 
11, 12, 13, 14, 15, 16, 17]. It has been noted earlier that pp interactions are only important 
for photospheric radii of GRB fireballs [18, 19]. The interactions of shock accelerated protons 
with low energy radiations (^7) lead to the production of ultrahigh energy mesons and leptons 
P7 —7- TT^'^X, vr"*" —7- fJ-^i^^, /i^ -^ e^v'^i'e, X can be neutrons and they will decay to protons, 
electrons and antineutrinos of electron flavour (n — > p + e~ + V^) in 886 sec in their rest frames. 
The other secondary products in p^ interactions are p'y — )• K~^'^X where X can be either A'^, T,^ or 
S"*". Kaons decay to lighter mesons, leptons and neutrinos. If the magnetic field inside the GRBs 
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is very high then the ultrahigh energy muons, pions and kaons lose energy significantly before 
decaying to neutrinos. Synchrotron cooling of charged muons and mesons would lead to decrease 
in the fluxes of neutrinos from their decay. Although, the cross-section for kaon production from 
P7 is much less than that for photo-pion production, kaons being heavier cool at a much slower 
rate compared to muons and pions. As a result at very high energy the kaon decay channel of 
neutrino production becomes significant compared to the pion decay channel [20, 21, 22, 23, 24]. 
Shock accelerated protons are expected to produce high energy neutrons in various interactions, 
and these unstable particles subsequently decay to leptons and protons. 

In the context of UHECR production from GRBs the neutron decay channel was considered earlier 
in [25]. The energy and time dependence of neutron, neutrino fluxes from GRBs were calculated 
within the scenario of external shock model in this paper. Also the radiation halos from neutron 
decay electrons, protons and their detectability in different wavelengths were studied. 
In a recent paper [26] the authors have explored the possiblity of detecting high energy diffuse 
neutrino flux from GRBs by IceCube, assuming the diffuse UHECR flux as detected by HiRes is 
the proton flux produced by the decay of cosmic ray neutrons in GRB fireballs. 
In our work we have not assumed that GRBs are the sources of the UHECR events observed 
by Pierre Auger [1, 2], HiRes [3] or AGASA [4] experiments. We have considered the possibility 
of cosmic ray acceleration in individual GRBs to ultra high energies to explore the parameter 
dependence of the high energy neutrino fluxes produced in decay of various particles. We discussed 
the importance of the neutron decay channel in our earlier paper [27] for cosmic accelerators with 
high internal magnetic fields. In this paper we have calculated the neutrino fluxes from GRBs in 
the internal shock model through muon, photo-pion, neutron and kaon decay, including the effect 
of proton energy loss by synchrotron cooling and P7 interactions, synchrotron cooling of charged 
muons, pions, and kaons in the internal magnetic field of GRB. After including proton energy 
losses, we find that the neutron decay channel of high energy antineutrino production may remain 
important in GRBs depending on the values of the GRB parameters. Moreover, the various decay 
processes may show distinct features in the neutrino flavour ratios. 

Due to the ongoing experimental activities to detect high energy neutrinos from GRBs [28, 29, 30, 
31, 32] this field has remained exciting. IceCube collaboration has claimed to reach the sensitivity 
of detecting neutrino flux from GRBs at TeV energy. IceCube opearted in a 40-string configuration 
from April 5, 2008 to May 20, 2009 [29]. They considered 117 bursts and no events were detected 
above the atmospheric background. For cosmic neutrinos with an i?"^ energy spectrum an integral 
flux limit of E'^cf) < 3.6 x 10^^ GeVcm^'^sec~^sr^^ has been found in the energy interval of 2 x 
10^ — 6.3 X lO^GeV [30]. In particular the non-detection of neutrinos from GRBs [32] has placed 
a tighter upper bound which is 3.7 times below the theoretical predictions [7, 8, 11, 26] combining 
40 and 59 strings. IceCube collaboration has concluded [32] that GRBs are not the only sources 
of cosmic rays above energy 1 EeV or the efficiency of neutrino production is much lower than the 
current predictions. The values of the bulk Lorentz factor T and the ratio of energies in protons to 
electrons have been constrained at the 90% confidence level in their Fig. 4. Hiimmer et al. [33] have 
recalculated the neutrino fluxes from GRBs and concluded that their result is significantly below 
the limit set by IceCube experiment operated in 40 strings. It has been pointed out by Li [34] 
that the theoretical prediction of neutrino flux from GRBs in IceCube papers is an overestimation, 
and the non-detection of GRB neutrinos is consistent with the correct theoretical estimation of 
neutrino flux. He et al. [35] have calculated the neutrino flux from GRBs including proton energy 
losses and they have also suggested that the neutrino flux predicted theoretically in the papers by 
IceCube collaborations is an overestimation. In calculating the photon number density IceCube 
collaboration approximated the energy of all photons by the break energy of the photon spectrum. 
The recalculated neutrino flux by He et al. [35] from 215 GRBs observed by IceCube 40 and 59 string 



configurations is about 36% of the 90% confidence upper limit obtained by IceCube collaboration 
and it is consistent with non-detection of GRB neutrinos by IceCube detector. The objective of the 
present work is to show that the relative importance of the different channels (pion, muon, neutron 
and kaon decay) of neutrino production in GRBs depends on the values of the GRB parameters. 
We have compared our calculated flavour ratios with earlier work [36, 23]. After including muon 
energy loss and kaon decay channels in detail at high energy, our flavour ratios are lower than the 
constant value of (pu^/i't'ue + ^P^t) — 0-64 found in earlier work. 

2 High energy neutrinos from gamma ray bursts 

Frames of references have been assigned as "c" for comoving or wind rest frame, "p" for proton 
rest frame. Quantities measured in the source rest frame are written without any subscript. In this 
paper we consider the muon, photo-pion, neutron and kaon decay channels of high energy neutrino 
production in the internal shock model of GRB prompt emission. We have used the low energy 
photon flux typically observed by Swift in the energy range of 1 KeV to 10 MeV to calculate the 
neutrino flux from individual GRBs. The photon energy spectrum from a GRB can be expressed 
as a broken power law with break at e^ in the source rest frame related to the break energy in the 
comoving frame e^ ,, as e^ 

J^ ( ^ -71 ^ ^ ^b 

(1) 

7i < 2, and 72 > 2. The normalization constant A is related to the internal energy density U by, 
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These photons are interacting with shock accelerated protons to produce charged and neutral pions 
and kaons. The maximum energy of the shock accelerated protons in the GRB fireball can be 
calculated by comparing the minimum of the p7 interaction time scale (^^7), p-synchrotron cooling 
time scale (tsyn) and dynamical time scale (tdyn) of a GRB with the acceleration time scale (tacc) 
of the protons as discussed in [37] . 

The inverse of the time scale of p^ interactions with protons of energy ep^c in the comoving/wind 
rest frame leading to the production of secondary particles 'a' is, 

taHe,,c) - -^^-^ (3) 
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where Tp^c = ^p,c/^rnpC^, and e^^th is the threshold energy of photon in proton rest frame, e^^p and 
e^^c are the photon energies in the proton rest frame and the comoving frame respectively. ia{^'y,p) 
is the average fraction of energy lost to the secondary particle by a proton of energy ep^c- The 



expression for the low energy photon flux can be substituted from equ.(l). The second integration 
gives, 
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72+1 2rp,c '=7.c 2rp,c ^ ^7,c 

The inverse of the time scale of p7 interactions leading to the production of secondary particles 'a' 
can be expresses as , 

^a^(ep,c) = 15T;2~ / '^e7,P'^(^7,p)?(^7,p)<^7,P-^2- (6) 

'^'' p,c •'<^-,,th 

The fractional energy transferred from proton to the secondary particle in an internal shock of 
radius r^ is fa = rd/{Tcta). 

2.1 Neutrinos from photo-pion decay 

We have considered production of pious in pj interactions through the decay of resonant particle 
A+. At the delta resonance both vr'' and 7r+ have been assumed to be produced with equal 
probabilities, vr^ gets on the average 20% of the proton's energy. The charged pious decay to 
muons and neutrinos. Finally the muons decay to electrons and neutrinos, antineutrinos. Each 
pion decay followed by muon decay gives two neutrinos, one antineutrino and one positron. If 
the final state leptons share the pion energy equally then each neutrino carries 5% of the initial 
proton's energy. The minimum energy of the protons interacting with photons of energy e?, for the 
A resonance with cross section ua -^ 5 x 10"^^ cm'^ in the source rest frame is, 

6^,A = 1-3 X 10^rioo(6!;,Mey)~'GeV. (7) 

The first break in tiie energy spectrum of neutrino, e^^ is due to the break in the photon spectrum 

at Cy. 

e' = 6.5 X W'rl,, (e^ ,,,,.) "' GeV. (8) 



The fireball Lorentz factor T^qq = r/300, photon luminosity ^7,51 = L^/ (10^^ ergs /sec), variability 
time tt,,-3 = (tt,/10^^sec) are the important parameters of a GRB. 

The minimum Lorentz boost factor of a GRB wind for observability of photon of energy et (energy 
measured in wind rest frame) depends on L^ and t^ as mentioned in [38]. 
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The internal energy density U relates to photon luminosity, L^ = Anr^'^T'^cU. r^ = T'^cty is the 
internal shock radius. The charged pions and muons lose energy due to synchrotron radiation in 
magnetic fields inside the sources before decaying to neutrinos. The expression for magnetic field 
in the comoving frame is, 

5e = 5.04 X 10^ ef_, l]% r^,l t^ig G (10) 

€b and Cg are the energy fractions carried by the magnetic field and the electrons respectively. Pion 
cooling energy is ten times higher than muon cooling energy. We would be overestimating the total 
neutrino flux if we use the pion cooling energy to derive the second break energy in the neutrino 



spectrum as the number of neutrinos produced from muon decay is more. The muon coohng break 
energy el „ can be expressed as a function of GRB parameters. 

el^^ = 2.56 X lQ^el/\-^^''L-\!M,,t,,_,GeY. (11) 

The total energy to be emitted by neutrinos of energy Cj^^Tr from photo-pion decay (considering 
muon and pion decay neutrinos together) in the source rest frame of a GRB is [13], 

,2 dN, ^ 3/. (1 - 6e - 6b) ^,so \ ^ _^ e^,- < <,M ..„. 

where El^° is the total isotropic energy of the emitted gamma-ray photons in the energy range of 
IkeV to lOMeV, which is available from observations. It is the product of L^ with the duration 
of the prompt emission from the GRB. In deriving eqn.(12) it has been assumed that the shock 
accelerated electrons in the internal shock model are cooling fast. They are in the radiatively 
efficient regime and their energy goes to the emitted photons. The fraction ep = 1 — eg — e^ of shock 
energy is assumed to be all in relativistic power law protons. This is an extreme assumption that 
the shocks are 100% efficient in converting bulk relativistic shock energy into relativistic power law 
electrons, protons and fields. However, it is known from simulations that in shocks many of the 
shocked protons remain in a thermal peak, and only some unspecified fraction of the total number 
of protons ends up in the power law. This reduces the power law proton (and neutrino) flux by a 
multiplicative proton injection fraction rjp < 1. 

The maximum energy of neutrinos from pion decay is approximately 5% of the maximum energy 
of protons (ep^max)- fn is the fractional energy of a proton going to pion inside the GRB fireball. 
It can be expressed in terms of the parameters of a GRB. 

Mep) = fS{ Y'-^L-i (13) 

1.34^2-1 / ep \ f ^ fb 

72+1 y^i^) ^p<^pa 

where f^ = C-kt^ — r — '^ t — j — ■ — j — r and ^^ = 0.2. k is a normalization factor in eqn.(12). The 

relativistic electrons produce the photons by synchrotron radiation and inverse compton scattering 
of low energy photons. Four orders of magnitude in photon energy corresponds to two orders of 
magnitude in the energy of the radiating charged leptons. The ultra-relativistic electrons have a 
power law spectrum, spectral index is assumed to be —2.5. This corresponds to k = 1.8 assuming 
photon fluence is proportional to neutrino luminosity. 

2.2 Neutrinos from neutron /3-Decay 

Ultrahigh energy neutrons are also produced in p^ interactions along with pions and kaons [25, 26]. 
This channel was considered to estimate TeV neutrino flux from Cygnus 0B2 located about 1.7 
kpc away from us [39]. We follow a similar formalism. The neutrons (with Lorentz factor r„) decay 
(n -^ p + e~' +T'e) to Pe with a decay mean free path cTnfn = W{en/EeV) Kpc. r„ = 886 seconds 
is the lifetime of a neutron in its rest frame and e„ is its energy in the source rest frame. The fraction 
of a proton's energy lost to neutron production in the process P7 — > iT~^n at the A— resonance [27] 
is, 

1.34^1-1 f JiL_y^ f ^ fb 
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where Jq = ^n-Fz — r — ^^4^^ t — j — - — j — r and ^„ = 0.8. The neutron spectrum can be expressed in 

terms of El^^° and the equipartition parameters eg, €b as follows, 

^2 dNn ^ fn l- ee-es ^iso ,^^. 

We have assumed the probability of production of neutrons in resonant p7 interactions to be half. 
The energy flux of antineutrinos of energy ep^n can be estimated with the neutron flux [dNn/ den) 
at the source rest frame [39] as, 



2 dNp 

^^ij n. 



^n.max 



den dNn ( Dsrrin \ m^ 



den \ / 2 eo 



X ep,n- (16) 



eo is the mean energy of an antineutrino in the neutron rest frame. The bracketed term is the 
decay probability for a neutron with energy e„ travelling a distance Dg. However if the source 
is at a large distance then all the neutrons are expected to decay before reaching earth. The 
maximum energy of the secondary neutrons is en,max = £,n^p,max- The maximum energy of the 
antineutrinos from neutron decay is Q/run times the maximum energy of the decaying neutrons, 
where Q = m„ — f^p — inrie and m^, rrip^ rUe are the masses of neutron, proton and electron 
respectively. 

2.3 Neutrinos from kaon decay 

In p7 interactions charged kaons {K^) are produced through the following interactions p7 — ?> 
K^ hS* and p7 -^ K^TP . Although the cross-sections of these interactions are lower compared to 
photo-pion production [40], at very high energy the total neutrino flux produced through kaon 
decay becomes higher than that from photo-pion decay [21, 22]. The cross-section for K^ and 
A'' production peaks at a\o x+ = 2 x lO"^'^ cm? for photon energy eg = 1-3 GeV in the proton 
rest frame and it has a full width at half maximum of 6ej^+= 0.9 GeV. K"^ and T,^ are produced 
with the highest cross-section of o"so^if+ = 2.4 x lO"^'^ crn^ at e^ = 1.45 GeV, and 6ex+ = 
0.8 GeV in this case. The average fraction of energy lost by a proton in producing a K~^, is 

-I m| (A0ori;0)-m2 

^K+ — 21^ s ^^\- S i^ t^^ invariance of the square of the total four-momentum of 

the p7 system and m\, m^, mj^+ are the masses of A*^, Ti^ and K'^ respectively. The threshold 
energy of photon in proton rest frame for kaon and A*^ production is e^/j x+ '^ 455 MeV. Fig. 1 shows 
the inelasticity for photo-pion and K~^ with A° production in p7 interactions plotted against the 
photon energy in the proton rest frame. The cross-sections are measured upto e-y^p = 2.6GeV in 
the paper by [40]. The fractional energy transferred from shock accelerated protons to kaons has 
been calculated in the same way as for photo-pions. 

In the case of resonant production of K~^ with A^, the fractional energy loss of a proton of energy 
ep is the ratio of kaon production and dynamical time scale of the GRB, 
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The protons of energy e^ . satisfy the threshold energy condition for production of K^ and A*^ in 



In the above equation Jq = .^x+Fi — 7 ''^ — t — i — ^ — i — T '^ith the mean value of ^x+ ~ 0-2 [21]. 
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Figure 1: Inelasticity of meson production from proton photon interaction as a function of photon 
energy in proton rest frame. 

interactions with low energy photons of energy e^^^gy. The expression for e^^o is as given below 
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If we consider p^y interactions with multiparticles at the final states, p^ -^ K'^A^tt^, then the cross- 
section for K^ production is cJx+ mult ^ 0-5 x 10"^" cm^ as measured by [41] for photon energy 
upto 2.6GeV in proton rest frame. In this case we have assumed 20% of a proton's energy goes to 
K+. 
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4,muit = 4.9 X 10^rioo(6^,Mey)~' GcV. (20) 

K^ decay to secondary neutrinos by the following channels, K~^ -^ ^"'"f^(63%), 7r^7r'^(21%), 
7r+7r+7r-(6%), 7r°e+i/e(5%), 7r^fj,+Uf,{3%) and 7r+7r07r0(2%). Due to their heavier mass K+ cools 
at higher energy compared to muon and pion. The synchrotron cooling break energy in the kaon 
spectrum is at 

e^ = 2.2 X W^el/%'^'L-l,{'Ti,,U,.s GeV. (21) 

It is derived by comparing the decay and cooling time scales of kaons. Although the cross-section 
of kaon production is much less than that of pion production, the neutrino flux from kaon channel 
exceeds the flux from pion channel at very high energy due to the slower rate of cooling of kaons. 



Similarly one can do the calculation for K'^ and TP production in p7 interactions. The neutrino 
flux from this channel will be added to the previous channel (eqn. 22). 

Neutral kaons are produced in p^ interactions with a cross-section a^o s ~ 0.6 x 10^^*^ cm?' [41] at 
the peak energy eo^ = 1.45 GeV and width de^o = 0.7 GeV. Half of the neutral kaons are assumed 
to be long lived kaons {K^). K^ can be produced in p^ interactions with multiparticle final states 
{Pl ~^ i^^A''7r+, i^£A*^7r+ and K'^Ti^tt^). The cross-sections of these interactions are measured 
as 0.5 X 10~^°cm^, 0.5 x lO"^*' cm^ and 0.2 x 10"^'^ cm^ respectively [41]. K^ decays through the 
following channels vr+e" 27^(39%), 7r+ n' i7f,{27%) , 7r°7r°7rO(21%), and 7r+7r-7r0(13%). K^ decays to 
two charged pions through this channel -7r"'"7r~(69%). The pions finally decay to neutrinos and 
antineutrinos. The total neutrino flux from K^ and Kg can be calculated in the same way as K^. 
The energy flux of neutrinos with energy eu,K from kaon decay in GRB source rest frame is, 






fu^K is the fractional energy of a kaon transferred to each neutrino/antineutrino produced in kaon 
decay. The observed total neutrino flux with energy e° on earth is, 



where z is the redshift of the GRB. 

2.4 Results and Discussions 

Our calculations are based on the standard internal shock model of GRBs. In internal shocks the 
shock radius r^ is related to the bulk Lorentz factor F and variability time t^, r^ = T'^cty, where c is 
the speed of light. We have not assumed any relation among the GRB parameters T and isotropic 
energy [42, 43, 44] or peak luminosity and observed break energy in the low energy photon spectrum 
[44]. The values of the parameters of GRBs have been varied to show that the importance of the 
different channels may vary depending on the luminosities, Lorentz factors, variability times and 
other parameters of GRBs. The redshifts of the GRBs have been denoted by z. We have assumed 
the GRBs have durations of prompt emission 5 sec. Our results for various set of values of the 
parameters of the GRBs are shown in Figure 2 to Figure 6. The following values of the parameters 
have been used. 

Figure 2. 71 = 1, 72 = 2.2, z = 1.8, Ly = 10^^ erg/sec, F = 600, t^ = 20 msec, e^ = 0.5Mey, 
eB=ee = 0.3, /° = 0.09 and r^ = 2.16 x 10^^ cm, B^ = 1.54 x 10^ G. 

Figure 3. 71 = 1.2, 72 = 2.5, z = 1.2, L^ = 10^^ grg/sec, F = 600, t„ = 20 msec, e^ = 0.5Mey, 
CB = 0.6, ee = 0.06, /° = 0.14 and r^ = 2.16 x 10^^ cm, B^ = 1.219 x 10^ G. 

Figure 4. 71 = 1.8, 72 = 2.01, z = 1, L^ = 5 x 10^^ erg/sec, F = 130, t^ = 25 msec, e\ = 0.5MeV, 
eB=ee = 0.3, /° = 0.16 and r^ = 1.26 x lO^^ ^m, B^ = 0.222 x 10^ G. 

Figure 5. 71 = 1.2, 72 = 2.2, z = 0.8, L^ = 10^^ erg/sec, F = 1000, U = 10 msec, e^ = 0.5Mey, 
eB=ee = 0.1, /° = 0.054 and r^ = 3 x 10^'' cm, B^ = 4.303 x 10^ G. 



Figure 6. 71 = 1, 72 = 2.2, z = 1, L^ = 10^^ erg/sec, T = 600, t„ = 2 msec, e^ = 0.5MeV, 
eB=ee = 0.3, /o = 0.009 and r^ = 2.16 x 10^'' cm, B^ = 1.543 x 10^ G. 

The total neutrino fluxes from the decay of muons, pions and neutrons are shown in Figure "a". 
The total neutrino flux from the decay of muons and pions, where the pions were produced in 
resonant interactions of shock accelerated protons with low energy photons, is shown by solid blue 
line. Red solid line is for the total neutrino flux from the decay of neutrons, where the neutrons were 
produced in resonant pj interactions. At low energy the total neutrino flux from neutron decay 
can exceed the sum of the total neutrino fluxes from muon and pion decay as shown in Figure 3 [a] 
and 3[c]. At very high energy the neutron decay channel may also become important as shown in 
Figure 4 [a] and 4[c]. 

The total neutrino fluxes from the decay of charged and neutral kaons are shown in Figure "b" . 
Charged kaons (K^) are produced in resonant interactions of shock accelerated protons with the 
low energy photons. The total neutrino fluxes produced in the decay of charged kaons are shown 
by solid orange line for the resonant production of K^. The dashed orange line is for the total 
neutrino flux from the decay of K~^ where the charged kaons were produced in pj interactions 
with multiparticle final states. Similarly solid, dashed cyan and magenta lines represent the total 
neutrino fluxes from decaying long lived (K^) and short lived (Kg) neutral kaons for the above 
mentioned cases. Green solid line shows the total neutrino flux from the charged and neutral kaons 
after summing up all the different components mentioned above. The total neutrino flux produced 
from kaon decay shows distinct peak at the highest energy. 

Figure "c" represents the total fluxes of neutrinos from the decay of different secondary particles 
(/x+, TT"'" (blue), neutron (red) and kaon (green)). After adding up the contributions from all the 
decaying secondaries the total neutrino flux has been shown by black solid line. 
The peak in the neutrino spectrum near PeV energy is due to protons of energy 200 PeV. 
The flavour ratio of neutrinos can be an useful observable to reveal the underlying physical phe- 
nomena in many processes [36, 45, 46]. Kashti & Waxman [36] explored the energy dependence of 
the neutrino flavour ratios produced in pp and p'j interactions. Energetic pions and muons cool by 
emitting radiations in the internal magnetic field of the sources. As muons cool faster compared to 
pions muon decay neutrino flux is suppressed at high energy. This leads to a flux ratio of neutrinos 
and antineutrinos (p^^ : (p,^^ : (p^^ = 1 : 1.8 : 1.8 expected on earth instead of 1 : 1 : 1. This flux 
ratio corresponds to R = 0^/((/>e + 0^) = 0.64 on earth. 

In the context of high energy neutrinos from GRBs the flavour ratio R = (/)^/{(/)e + (pr) expected on 
earth has been calculated in [23]. In our Figure "d" we have plotted the flavour ratio R expected 
on earth from our calculations. Our plots can be compared with those obtained by [23]. Due to the 
contributions from neutron and kaon decay neutrinos and antineutrinos the plots of flavour ratio 
R in our work show different features. In Figure 4[d] the sudden dips are due to the importance 
of the neutron decay channel (near 10 GeV) and kaon decay channels (near 10^ GeV) compared 
to the pion decay channel. In rest of the plots (Figure 2[d], 3[d], 5[d], 6[d]) showing the energy 
dependence of R, the first bump is due to the decrease in neutron decay antineutrinos above 10^ GeV 
and the second bump is due to muon cooling. The step at higher energy after the second bump 
is due to the importance of the kaon decay channel compared to the pion decay channel. In the 
paper by Baerwald et al. [23] below 100 TeV the flavour ratio R is between 0.36 and 0.5 after 
including variations in values of all the GRB parameters. Above 10 PeV they have R = 0.64 due 
to muon damping assuming at the source the flux ratios are (p^^ : (p^^j^ '■ <pu-r =0:1:0. We have 
not made such assumptions in our calculations. We have included muon cooling at high energy 
with the neutrino flux given in eqn.(12). We have also included the kaon decay neutrinos in our 
calculation of R, as a result our calculated R does not remain constant at high energy. We have 



used sin^20i3 = 0.1, where ^13 is the neutrino mixing angle measured by DOUBLE CHOOZ [47], 
DAYA-BAY [48] and RENO collaborations [49]. In our Figure 2[d] at 5 TeV the flavour ratio is 
R = 0.41, however in Figure 3[d] neutron decay channel is more important at low energy and 
R = 0.36. These results can be compared with a pure neutron beam source where R = 0.28 and for 
a pure pion beam source this ratio is 0.5. In Figure 4[d] upto 1 PeV the flavour ratio R is 0.47 then 
it goes down to 0.37 at 7.5 PeV as the neutron decay channel dominates at this energy. In Figure 
5[d] and 6[d] at 1 TeV the flavour ratio is R = 0.41. In Figure 2[d], 3[d], 5[d] and 6[d] near 1 PeV 
muon and pion decay neutrino flux becomes higher than the neutron decay antineutrino flux and 
R increases to 0.48. In Figure 2[d] due to muon cooling R increases to 0.54 near 2 EeV and then it 
attains the highest value of 0.57 at 10 EeV due to the kaon decay channels. In Figure 3[d] near 1 
EeV muon cooling starts, at the same time kaon channel becomes important. R is 0.51 at 1 EeV 
and near 30 EeV R is 0.57. In Figure 4[d], 5[d] and 6[d] this highest value of R is attained above 1 
EeV, 80 EeV and 10 EeV respectively. At 1 EeV energy our R values are almost 20% to 30% lower 
than the value of 0.64 obtained by Baerwald et al. [23]. At the highest energy our R values of 0.57 
are 12% lower than the constant value of 0.64 in Baerwald et al. [23]. One may notice from our 
figures that the muon cooling and the importance of the kaon decay channels together determine 
the ratio R at high energy. 

It is interesting to note that the flavour ratio R carries the information of the origin of the neutrinos 
and the antineutrinos. 

We have shown that the different decay processes may contribute differently depending on the 
values of the parameters of the GRBs. The energy of the neutrinos above which the neutron decay 
channel will dominate over the pion decay channel can be derived by equating the fluxes from the 
two channels with eu,n=Gu,n- We compare the total neutrino spectrum from pion decay above the 
synchrotron cooling break energy (eqn. 12) with the same from neutron decay (eqn. 16) assuming 
all the neutrons have decayed during propagation. 



eu' 






'mneu,n/'^^G 



> 1.96x10^^(3-71) ( ^ 1 ( ^^^ 1 ( ^ {'"' et,,,,/GeV 




(24) 

The energies 6-?g.33 = erViW^'-'^'GeV), 6^,^,6.4 = etJ{W^-''GeV), e^j,,, = ei/{W'-^' GeV) 
and e[^^5gi = e^^/(10^'®^ GeV) are calculated for the following set of GRB parameters 71 = 1.2, 
72 = 2.3, Lj = 10^^ erg/sec, F = 300, t^ = 1. msec, e'^ = IMeV. Solving this equation numerically 
one can find the energy at which neutrinos from pion decay channels will be suppressed below the 
neutron decay channel. Similar to eqn. (24) one can get the inequality relation for the energy 
above which the neutrino flux from the kaon decay channel (eqn. 22) will dominate over the pion 
decay channel (eqn. 12) 

e.^ > 4.09 X 10^^ (1.54)-- V.X f-l^') f^^^V' '<,,6./Gey^ (25) 



,^0,0.001/ \ ''p,A,7.11 



b 
p,A0,7.59 



where f,y^k is the fractional energy of a kaon going to each neutrino/antineutrino, and e[ 
el^Ao/iW^-^^GeV), ep,A,7.ii = ep,A/(10'^'" G'el/), same set of values of the GRB parameters has 
been used here. These inequalities depend on the GRB parameters, L^, F, t„ and 71, 72, break 
energy in the photon spectrum. A small change in values of any of the parameters will lower one 
channel with respect to another and change the shape of the total neutrino spectrum from a GRB. 
The effect of stochastic reacceleration on photon and neutrino spectra from GRBs has been studied 

10 



for photospheric radii of GRB fireballs and high Lorentz boost factors in [50] . Inclusion of stochastic 
reacceleration may shift the peak in the total neutrino fluences to a higher energy as given in their 
Fig. 7. In that case the contribution of the neutron decay channel to the total neutrino spectrum 
may still remain important at low energy. We have discussed about the limits obtained by IceCube 
experiment in the Introduction of this paper. We compare our results with the recent IceCube limit 
[32] for GRB neutrinos with 40 and 59 strings. The details of their calculations and the values of 
the parameters they have used are given in [28, 29, 32]. IceCube collaboration observed 117 GRBs 
with 40-string, two additional GRBs included from test runs before the official start of 59-string 
and 181 GRBs with 59-string. The sum of the neutrino energy fluences from all the 353 GRBs 
observed is shown in their Figure. 1. of [32] with the lower most black solid line. We have shown 
in our Figure [c] the same with solid brown line [32]. Our calculated GRB neutrino fluences from 
353 GRBs for each set of parameter values are shown with brown dashed lines in Figure 2[c]-6[c]. 
Our choice of values of the GRB parameters give neutrino fluences below the current limit set by 
the IceCube expriment. 

3 Conclusion 

Gamma ray bursts have been speculated to produce high energy particles like cosmic rays, gamma 
rays and neutrinos. Ultrahigh energy muons, pions, neutrons and kaons could be produced in GRBs 
due to p7 interactions. These particles subsequently decay to produce very high energy neutrinos 
and antineutrinos. Charged particles lose their energy due to synchrotron cooling in the magnetic 
field before decaying to lighter particles. Kaons suffer synchrotron loss at higher energy compared 
to muons and pions. As a result at higher energy the kaon decay channel of neutrino production 
dominates over the muon and pion decay channels. We have included all the decay channels of 
charged {K~^) and neutral kaons {K^ and Kg) to study this effect. The total neutrino flux from 
neutron decay channel may be more important than that from muon and pion decay channels 
depending on the values of the parameters of GRBs. The flavour ratios of neutrinos show distinct 
features depending on the values of the GRB parameters and the dominance of different decay 
processes at different energies. IceCube has set limit on the GRB neutrino flux. Future detectors 
with higher sensitivities will be able to reveal the role of GRBs as a UHECR accelerator. 
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Figure 2: 71 = 1, 72 = 2.2, z = 1.8, L^ = lO^^erg/sec, T = 600, t^ = 20msec, e^ = 0.5MeV, 
e^ = Ce = 0.3. [a] Blue solid: resonant pion, muon, red solid: resonant neutron decay neutrinos; 
[b] orange solid, dashed: resonant, multiparticle final state K^ decay neutrinos, cyan solid, dashed 
for i^£, magenta sohd, dashed for K^, green sohd: total neutrinos from kaon decay; [c] blue: total 
neutrinos from muon and pion decay, red: neutron decay, green: kaon decay; black: total neutrinos 
from decay of all particles; brown solid: IceCube-40+59 string limit, brown dashed: our aggregate 
of 353 GRBs [d] flavour ratio R as discussed in the text 
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Figure 3: 71 = 1.2, 72 = 2.5, z = 1.2, L^ = lO^^erg/sec, T = 600, U = 20msec, e° = O.SMeV, 
eb = 0.6, ee = 0.06. Line styles same as in Figure 2. 
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Figure 4: 71 = 1.8, 72 = 2.01, z = 1, L^ = 5 x lO^^erg/sec, T = 130, t„ = 25msec, e° = O.SMeV, 
€5 = ee = 0.3. Line styles same as in Figure 2 



14 



10" 



10'' 



10- 



~i 1 1 1 1 1 1 r 



[a] 




_] I I I L 



£?* [GeV] 





€!'lGeV] 



0.7 


1 1 1 1 1 


1 1 

[d] 


1 


0.6 


- 




- 


0.5 




J 




0.4 


1 1 1 1 1 


1 1 





<* [GeV] 



€!'lGeV] 



Figure 5: 71 = 1.2, 72 = 2.2, z = 0.8, L^ = lO^^erg/sec, T = 1000, t^ = 10msec, e^ = O.SMeV, 
cb = Ce = 0.1. Line styles same as in Figure 2. 
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Figure 6: 71 = 1, 72 = 2.2, z = I, L^ = lO^^erg/sec, T = 600, t^ = 2msec, e\ = O.SMeV, 
es = Ce = 0.3. Line styles same as in Figure 2. 
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